Correlations of superconducting fluctuations in a two-gap system 
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We derive analytically the spatial correlation functions for gap fluctuations in two-band scenario 
with intra- and interband pair-transfer interactions. These functions demonstrate the changes in 
spatial functionality due to the presence of two channels of coherency described by the divergent 
and finite correlation lengths. Even at the phase transition point both channels are essential for 
two-band superconductivity. Generally their relative contributions depend on the temperature and 
system parameters. 
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I. INTRODUCTION 

The theory of superconductivity with overlapping 
bands has started to develop since 1959^ however, only 
after discovery of multicomponcnt nature of MgB2 in 
2001— and pnictides in 2008 3 the multigap approaches 
have become an object of exceeding interest. 

The peculiarities of the spatial coherency in multiband 
superconductors have attracted much attention recently 
in connection with type- 1.5 behaviour—. In usual one- 
band systems there is only one coherence length which 
value in units of penetration depth determines response 
to magnetic fields, type-I or type-II. It was suggested 
that in two-band superconductor one has two correlation 
lengths resulting much richer physics than type-I/type-II 
dichotomy. In particular, there is a possibility to observe 
a mixture of domains of Meissner state and vortex clus- 
ters, called type-1.5 superconductivity. The latter regime 
is supported also by interband proximity effect and by 
different kinds of intercomponent interaction involving 
Josephson, mixed gradient or density-density couplings 6 . 

The existence of two qualitatively different length 
scales in two-band system was demonstrated more than 
twenty years ago 7 - and recently-^. Two distinct corre- 
lation lengths are also present in negative-iJ Hubbard 
model of two-orbital superconductivity-i£. In this respect 
the connection between peculiarities of spatial coherency 
and excitation of the Leggett mode in two-gap material 
was discussecUi. 

Different point of view on the correlation behaviour in 
two-band model is based on the statement that within the 
Ginzburg-Landau domain both order parameters vary on 
the same length scal e 12 ' 13 . An extension of the temper- 
ature domain indicates that two gaps are generally not 
proportional to each other, thus their spatial scales are 
decoupled only away from critical points. Numeric esti- 
mations for the healing lengths of the gaps confirm that 
conclusion for several superconducting materials^. Here 
we note that scientific discussion about discrepancy of 
length scales in the vicinity of critical temperature still 
continues^ 6 -. 

Experimentally, two characteristic length scales in var- 
ious two-band compounds were evidenced e.g. by vortex 



imaging^ 7 -, muon spin relaxation measurements^, in heat 
transport features^ as a function of magnetic field. 

In the present contribution we derive correlation func- 
tions for gap fluctuations in the two-band scenario. The 
spatial behaviour of these characteristics reveal two dif- 
ferent correlation lengths describing joint superconduct- 
ing condensate as a whole. These length scales are anal- 
ysed as the functions of the temperature and interband 
interaction constant. The competition between the con- 
tributions of the corresponding coherency channels to the 
correlation functions are discussed. 



II. DERIVATION OF CORRELATION 
FUNCTIONS 

We start with two-band superconductivity Hamilto- 
nian 

H = e a (k)a+ ks a Qk;s 

aks 

~V ^ ^ ^ aQ ' a akt a a-k+qi a a'-k'+qia«'k't> (1) 
aa' kk'q 

where e Q = e a — ji is the electron energy in the band 
a = 1, 2 relative to the chemical potential /i; V is the vol- 
ume of superconductor and W aa > is the matrix elements 
of intraband (a = a') or interband (a ^ a') pair transfer 
interaction. It is supposed that the chemical potential 
is located in the region of the bands overlapping. We 
assume that (effective) electron-electron interactions are 
nonzero only in the layer [i ± fiUD and W aa > is indepen- 
dent on electron wave vector in this layer. For simplicity 
we take W\% = W21. 

We calculate the partition function Z = Spexp {j^j^j 
for the macroscopic system by means of Hubbard- 
Stratonovich transformation 2 ^. For W 2 = W11W22 — 
W12 > and for real order parameters S a the static path 
approximation reads as 
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Here integration variables are treated as real and imagi- 
nary parts of Fourier components for non-equilibrium or- 
der parameters <5 Q (r) = ^ k S a ke lkr , F is non-equilibrium 
free energy of inhomogeneous system and F — F n in the 
absence of superconductivity. We do not expand the co- 
efficients 
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T c , which allows us to apply free energy in the 



and K n 
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form ^ substantially farther from critical temperature 
T c . Note that the coefficient K a used alludes isotropic 
situation. 

The homogeneous equilibrium state is defined by the 
minimization jg-\ s _ A = 0, which gives us the set 
of equations for coupled homogeneous order parameters 
Ai : 2, namely 



a a A a + b a A 3 a = cA 3 



(5) 



One should also take into account the relation between 
phases in equilibrium sqn(cAiA2) = +1. The criti- 
cal point is determined by condition ai(T c )a,2(T c ) = c 2 , 
which has two solutions T c ± and T c _ > T c+ . If T ca 
are intrinsic transition temperatures in the bands and 
Td > T c2 , then for Wi 2 -> we have T c _ -> T cl and 
T c+ — > T c2 . Note that in the system with coupled bands 
there is only one phase transition point T c = T c _. 

Now we linearise functional F near homogeneous state 
with free energy F^ by assuming S a (r) = A a + r\ a (r) . By 
means of complex Fourier components ry^k we have 



k>0 

where — A a + K a k 2 and A a = a a + 3b a A 2 t > 0. 
Note that due to interband pairing there appear non- 
diagonal terms in the quadratic form (jSJ). Statistics for 
the equilibrium state fluctuations is determined by the 

F 

distribution function e fc B T normalized to Z. By us- 
ing Gaussian approximation (J6j) we calculate mean values 
(Votk.'ila'ic) an d then correlation functions r QQ '(r — r') = 

Ek( f ' ak '(a'k) e!k ' r " r ' f° r the order parameter fluctua- 
tions considered at different points separated by distance 
|r - r'| ^ 0. We obtain T aa , = T+ a , + T~ a ,, where 
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Note also that r i2 = r 2 i- In the Eqs. ©-© we have 
introduced ^ = ^ and the correlation lengths £± are 
given by 



e ± 



mi 



e 2 ± 
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These quantities have the following properties. For finite 
interband pairing £_ > £ + > 0. In the temperature re- 
gion where £i > £ 2 one has £_ > £i and £ + < £ 2 . For 
opposite case £i < £ 2 we get £_ > £ 2 and £+ < £i. As 
a result, r^ a > 0. However, sqn(rf 2 ) = =Fsqn(c), i.e. 
depending on the sign of interband constant one contri- 
bution in Ti 2 becomes negative. 

The characteristics £± define the size of the region, 
where the order parameter fluctuations are significantly 
correlated. In fact, these length scales appear in the ex- 
ponents despite the band index taken for the correlation 
functions, i.e. £± describe joint superconducting state 
rather than individual bands. We note also that £± co- 
incide with the correlation lengths- found by means of 
inhomogeneous gap equations. 



III. RESULTS AND DISCUSSIONS 
A. Correlation lengths 

The presence of interacting order parameters makes 
the coherence properties of the two-band system quite 
different from the corresponding characteristics in single- 
band superconductors. To analyse the physics of one- 
band case one should take c — > 0. In this limit £± — > 
£a|c=0i i- e - one obtains two separate correlation length 
attributed to the band a = 1,2. Each length diverges 
at its own point given by intrinsic transition temper- 
ature T ca . Note that £_ — » £i| c =o and — > £ 2 | c =o 
in the temperature region where £i| c =o > £2|c=0j how- 
ever, £_ — > £ 2 |c=o a nd — » £i| c =o for the temperatures 
where £i| c =o < ^2|c=o- Further we assume for specificity 
T c2 < T c i, i- e the condition £i| c =o < ^2|c=o corresponds 
to the lower temperature region, while £i| c =o > ^2|c=o to 
the higher temperatures in the superconducting state. 

Non-zero coupling between bands modifies drastically 
the trivial physics of two non-interacting condensates. 
The coherency is described by lengths £± which become 
tricky combinations of band characteristics £ a , see Eq. 
©. To illustrate the evolution of £± with model pa- 
rameters we fix intraband ones: W\\ : 22 = 0.3 eVcell, 
pi, 2 = (1,0.94) (eVceirr 1 , u F1 , 2 = (5,5.104) • 10 5 m/s, 
ceil = 0.1 nm 3 . For these values T c2 = 0.81T c i. We also 



assume parabolic electron spectrum where ^ 
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FIG. 1. Above: The plots of the gaps A a as a solution of Eqs. (JS| (solid) and derived microscopically (dotted) vs temperature 
for various interband couplings Wi2- Below: The log plots of £_ (solid red), (solid blue), l\ (dashed red), £2 (dashed blue), 
wi(70%) (dotted red) and W2(70%) (dotted blue) vs temperature for same W\%. 



Fig. [T] shows temperature dependencies for correlation 
lengths together with the evolution of homogeneous gaps 
calculated numerically as interband coupling increases. 
We see that £_ and £+ as functions of the temperature are 
remarkably different. First, the length £_ behaves criti- 
cally diverging at phase transition point T c . At the same 
time £ + remains finite. Second, £_ can change below T c 
very non-monotonically, while the temperature depen- 
dence of £+ is substantially weaker^. The appearance 
of additional maximum in superconducting phase for £_ 
is strongly supported by the smaller values of W12, rep- 
resenting the memory effect about criticality in the band 
a = 2. The position of this maximum is correlated with 
the inflection point of the smaller gap which takes place 
in the vicinity of T C 2- As was pointed earlier—, the non- 
monotonicity of the critical coherence length elucidates 
the temperature behaviour of the gaps healing length 22 
and vortex size 2 ^ in a superconductor with weakly inter- 
acting bands. 

One can argue that the scheme based on the expansion 
© is applicable only close to critical point. We note that 
the coefficients (jj) taken allow us to go essentially farther 
below T c . For the comparison we have plotted in Fig. Q] 
homogeneous gaps calculated numerically by means of 
microscopic theory. The latter are approximated by the 
solutions of system ([5]) very well in the temperature re- 
gion considered. 

Due to the definition of the critical point 
ai(T c )a2(T c ) — c 2 and the relation A a (T c ) = a a (T c ) one 
obtains 



g?(rc)g 2 2 (T c ) 
£?(Tc) + ^(t c )' 



(10) 



and zero for the denominator of £_(T C ), i.e. the latter 
length diverges precisely at T c . This implies that only 
length scale £_ can be attributed directly to the super- 
conducting phase transition in a two-band model. In the 



vicinity of T c we get 
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one meets also in the 



k 2 pia 2 (T c )+p 2 ai(T c ) T c 

The expressions ([10 ]> -([TT j) 
literatur o 7 ' 12 . 

Next we denote the factor in Eq. (fTTj) by £ 2 (0), 
the value of the formula (jlll) at T = 0, and analyse 
£_(0) and £+(T c ) as the functions of interband inter- 
action. Fig. [5] shows these dependencies for different 
sets of intraband parameters. Analytic consideration in- 
dicates that £_|_(T C ) always decreases with | Wi 2 1 , while 
£_ (0) can pass through a maximum at some finite value 
of W±2- We interpret this feature as follows. The one- 
band limit c = gives T c = T c \ and a\{T c \) = 0. As 
a result, £ 2 | c=0 = Kl ^ l] T ^ T for T < T cl . This 
is standard one-band expression for the squared corre- 
lation length expanded near critical point. The factor 



e(o)i c =o 



ifi(T cl 



is proportional to i.e. £_(0)| c= o 



2pi j^v^v^^v^,* T £ 
decreases with the critical temperature increase and vice 
versa. In two-band system T c always grows with W12 
(see Fig. [5]) and one can expect the reduction of £-(0) 
with increase of \Wi2\ by analogy with one-component 
case. However, in two-component situation, especially 
for weak interband couplings, the memory effect related 
to the lower intrinsic phase transition is strong. The lat- 
ter is characterized by the temperature T c+ which always 
decreases with | W12 ] (see Fig. By analogy with one- 
band case it can lead to the rise of (0). Thus, there are 
two opposite tendencies associated with the temperatures 
T c ± which govern the behaviour of £_ (0) as a function of 
interband coupling. By analysing this competition ana- 
lytically we find that, if 
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FIG. 2. Left: The plots of £_(0) (solid) and £+(T c ) (dashed) 
normalized to their values at W12 — us interband coupling 
Wia. Right: The plots of T c = T c _ (solid) and T c+ (dotted) vs 
Wi2- Red curves correspond to the parameters given in text, 
blue ones to the modified parameters in band a = 2, namely, 
W22 = 2.26 eVcell, p 2 = 0.13 (eVcell)" 1 , v F2 = 10 6 m/s. 
For these values we have the same ratio T C 2 = 0.81T c i. 



£_(0) has maximum, whereas for opposite sign in Eq. 
(fT2| the function £_(0) has minimum at W\% = 0. We 
believe that non-monotonicity of (0) is clear footprint 
of the two-band nature near critical point. 

One comment should be made about non-critical co- 
herence length . The quantity is always finite and de- 
creases as the strength of interband interaction increases, 
crossing zero at W 2 = 0. At the same time, there is 
natural lower bound for coherence lengths in Ginzburg- 
Landau theory defined by the microscopic length scales 
(Cooper pair size in the bands) . The latter guar- 



antees the smallness of the gradient term in Ginzburg- 
Landau expansion. To estimate the maximal value of 
we use Eq. (fTtlf for c — 0. We find 



£+(T c )\c=o — £2pcl)|c=0 



1 



y/piWu - p 2 W : 



22 



(13) 



Consequently, the value of £+(T c ) can be magnified when 
T C 2 approaches T c \. In this process non-critical coher- 
ence length can surpass microscopic lengths^, i.e. two 
length scales of coherency found are meaningful even in 
the standard two-band Ginzburg-Landau model for rele- 
vant parameters. To overcome the restriction related to 
the microscopic lengths one should take into account the 
higher terms of the gradient expansion in the Ginzburg- 
Landau approach. In this way one gets better agreement 
with microscopic theory-2i. However, the theory based on 
two-band Eilenberger equations also predicts the disap- 
pearance of non-critical length for strong interband pair- 
ings at W 2 ~ OA The absence of the real non-critical 
correlation length may signal spatial periodicity of fluc- 
tuations of two-gap superconductivity^. 



B. Correlation functions 

Interaction between bands results more complicated 
structure of correlation functions as compared to the case 
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FIG. 3. The log plots of T aa i (dots) together with correspond- 
ing contributions F~ a , (red), F^a and |rf 2 | (blue) vs distance 
r - r' for T = 0.86T C (left column) and T = T c (right col- 
umn). Here W12 = 0.01 eVcell and intraband parameters as 
in text. At T c we have |±- « 2.2. The regimes "sd" and "Id" 
are discussed in text. 



of decoupled bands for which 



knT 



87r#Jr-r'| 



|r-r'| 
- e €<»lc=0. 



0. (14) 



Next we discuss the correlation functions for non- 
vanishing interband pairings. 

First, we consider different spatial regions. For shorter 
distances |r — r'| <gc ^ + < ^_ (denote as "sd") we obtain 
from Eqs. ©-© 



r 



sd 
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In this fully correlated case the functions T aa i are max- 
imal. If £i(T c ) 3> &(T C ), then near T c the main contri- 
bution to Tfl stems from critical, while to from non- 
critical channel of coherency and vice versa. Note that 
condition £i(T c ) 3> &{T C ) is supported by the smaller 
interband interaction. 

For larger distances £ + < <C |r — r'| (denote as 
"Id") the functions r^ Q , are defined mostly by critical 
contributions and they vanish. By approaching T c we 
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have in this regime £+ <C | r — r' | <C £- an d 
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Thus, at critical point Ti2 changes in space from constant 
value to the function which decreases linearly 

with log ^ r ' . The disagreement between K, a and K a 
characterises the behaviour of r QQ . If £,\{T C ) 3> £, 2 (T C ), 
then and T 1 ^ are the very same function at T c , but 
there is remarkable difference in the dependencies r 2 2 
and r 2 2 related to the change of the dominant coherency 
channel from non-critical to critical one. This transfor- 
mation is also noticeable in Fig. [31 and it is supported 
by the weaker interband couplings. For opposite situa- 
tion £i(T c ) <C £2(^0 ) we have different dependencies for 
r 8 ^ and rft, but same for Tf 2 and T 2 d 2 . Therefore, the 
changes in spatial functionality of the correlation func- 
tions are intrinsic for two-gap superconductors. 

To estimate the efficiency of different correlation chan- 
nels we find the crossover points where r+ a = T~ a and 
crossover width. For Tn the crossover takes place at 
distance 



is vanishing, i.e. non-critical channel dominates for all 
reasonable distances. 

r+ 



Finally, the relative contribution 



decreases and 



increases with distance. If £ a > these functions 



cross at £ n 



i.e. 5^ > tt 2 ^ for tiny |r — r' 



We define 



r QQ — r Q „ 

the width w a of crossover region as the size of the spatial 
area around £ a where and ffi^ simultaneously do 
not exceed fixed percentage p > 50%. We find 



w a {p) 



£+£- 



In- 



(20) 



At T c one obtains w a ~ i.e. the crossover between 
contributions in T aa takes place on the distance defined 
by non-critical length scale. The width w a shrinks at T c 
with increase of interband coupling. For nearly decou- 
pled bands wi ~ holds also in the vicinity of T c2 (see 
Fig- HI)! however, w\ grows with W12 at that tempera- 
tures. In Fig. [T]we distinguish clear crossover from close 
coexistence of contributions up to £&. The latter occurs 
in the region around Ti which widens as interband inter- 
action increases. 



IV. CONCLUSIONS 
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where £ 2 is the corresponding parameter derived for r 2 2- 
For fixed temperature the crossover takes place only in 
the behaviour of certain correlation function, T\i or T22- 
From definition (|18l) follows that there is the tempera- 
ture Te defined by the condition £i(T£) — & (Ti) where i a 
goes to zero. The position of Tg, is sensitive to the model 
parameters. If ufi > WF2 one has Tg < T c . However, for 
v-pi < vf2 there is a value 
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for which Tg = T Cl and for stronger interband interaction 
T t > T c . Next, for £ X (T) < £ 2 (T) we have £1 > and 

vice versa. At T c one obtains £2 = f+ln vh- Fig. [T] shows 
also that l a can substantially exceed especially for 
nearly decoupled bands. In this case the value T aa (£ a ) 



The evolution of correlation functions for two-band 
superconductivity indicates the presence of two distinct 
channels of coherency described by the critical (divergent 
at critical point) and non-critical (finite at critical point) 
correlation lengths. Although these characteristics are 
not related directly to the bands, two-component nature 
manifests itself e.g. in the non-monotonicities of crit- 
ical length scale as a function of the temperature and 
the strength of interband interaction. The features of 
the competition between coherency channels involved de- 
pend on the temperature as well as model parameters, 
e.g. coupling between bands or Fermi velocities. 
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